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Abstract 

Employing the Worldline casting of the Dosch-Simonov Stochastic Vacuum Model 
(SVM) for QCD, a simulated meson-meson scattering problem is studied in the Regge 
kinematical regime. The process is modelled in terms of the 'helicoidaF Wilson contour 
first introduced by Janic and Peschanski in a related study based on an AdS / CFT-type 
approach. Using lattice supported estimations for the behavior of a two-point, field 
strength correlation function, as defined in the framework of the SVM, the reggeon 
slope and intercept are calculated in a semiclassical approximation. The resulting 
values are in good agreement with accepted phenomenological ones. Going beyond this 
approximation, the contribution resulting from boundary fluctuations of the Wilson 
loop contour is also estimated. 



1. Introductory Remarks 

In addition to confinement, which constitutes a profoundly non-perturbative, problem 
and whose solution is of quintessential importance, for fully establishing QCD as a funda- 
mental theory for the strong interaction, there do exist specific dynamical processes, whose 
theoretical confrontation also calls for non-perturbative methods of analysis. One such situa- 
tion arises in connection with the theoretical description of high energy scattering amplitudes 
for which the soft sector of the theory is involved. From the experimental point of view, 
one such case arises in connection with Regge kinematics, entering directly the theoretical 
description of, among others, diffractive and low-x physics processes. In this paper we shall 
apply, in this specific context, the Field Strength Correlator Method [1], in the framework 
of the S VM, as has been formulated in the preceding paper (I) , i. e. in terms of its Worldline 
casting. In particular, we shall study a simulated case of a meson-meson scattering process 
whose quark-based description is of the general form 

(11) + (22) -> (33) + (44) 

adopting a standard picture, already employed in the QCD literature -see, for example [2] and 
[3], according to which quark 1 from the first meson and antiquark 2 from the second meson 
are very heavy, in comparison to the incoming, total energy s -hence their worldlines are 
considered to remain intact from the gluon field action. In turn this means that they can be 
described in the framework of the eikonal approximation. The light pairs 1, 2 and 3, 4, on the 
other hand, are annihilated and produced in the t-channel, where the eikonal approximation 
is not valid and a full treatment is called for their description. In the Worldline framework 
the process is schematically pictured in space-time by the straight eikonal lines (1 — > 3) and 
(2 — > 4), describing an intact quark and anti-quark and by the curves (1 — > 2) and (3 — > 4) 
which correspond, respectively, to the annihilated and produced quark antiquark pairs. The 
structure of the field theoretical amplitude can be written as follows , see Fig. 

G(x^x z ,x 2) xi) = (iS F (x 4 ,x 3 | AiS F (x 3 ,x 1 \ A)iS F (x 1 ,x 2 \ A)iS F (x 2 ,x A \ A)) a- (1) 

In the above expression iS F is the full fermionic propagator which, in the framework of the 
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Worldline formalism, assumes the form [4] 



iSp(y,x | A) = / dLe 



—Lm 2 



Dx{r)e 



x(0)=x 
x(L)=y 



m — 



$ (1/2) (L,0)Pexp (i J drx-Aj , 



(2) 

where <3> <J, - ) is the so-called spin factor (see paper I) for the matter particles entering the 
system. For us, it means that j — \. 

Inserting the above formula into Eq. (1) we find 

4 oo 

G(x 4 , x 3 , x 2 , xi) = jlj dr t 6(r t - r^e'^-^ j J Dx(r)5[x(r 3 ) - x 3 ]5[x(r 2 ) - Xl ] x 



i=i 



x(T4,)=X4, 



x5[x(ti) — x 2 ] exp 



1 T4 

— - J drx 2 (r) (spin structure) ^Pexp j> dx ■ A 



A 



(3) 



where the term spin structure corresponds to the following expression 



(spin structure) = \\ 



TO* - ^7 • X ( T i) 



^(t^O, (r = 0). 



(4) 



In principle, the Wilson loop appearing in Eq. (3) incorporates the dynamics (perturba- 
tive, as well as non-perturbative) of the process. In the framework of the Stochastic Vacuum 
Model (SVM) it assumes the form (see I) 



P exp j> dx ■ A 



exp 



= e~ A W. 



(5) 



~ / dS^z) j dS Xp (z')A^ Xp (z - z') 

S(C) S(C) 

The task to be undertaken in the present paper is to calculate the amplitude (3), us- 
ing the above expression which, it is reminded, gives the structure of the Wilson loop in 
framework of the SVM. The particular strategy to be adopted in our relevant effort can be 
outlined as follows. In Section 2 we perform a semiclassical calculation based on a combined 
minimization of the action -A[C], see Eq. (17) of I, with respect to the surface S[C] and of 
the surface S[C] with respect to the boundary C. This procedure will allow us to determine 
the dominant contribution to the Worldline integral (3) in the stochastic limit T^y/A 1, 
as determined in I. 
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As also noted in I, the first order approximation of the action A[C], in terms of the 
correlation length T g , is essentially the Nambu-Goto string action. The next order corrections 
give rise to terms which reveal different geometric characteristics of the surface S[C], given 
its embedding in a 4-dimensional background, such as, e.g., its extrinsic curvature. The 
presence of these terms, the origin of which is completely different from the known quantum 
corrections of the Nambu-Goto string, points out the powerful structure of the SVM. 

In Section 2 we proceed further to take into account the rigidity of the surface S[C], 
which, as it will turn out, plays an important role for determining the reggeon intercept. Of 
lesser importance, but in any case computable, are the corrections related to the fluctua- 
tions of the (Wilson) 'curve' which forms the boundary of the surface and will be discussed 
in Section 3. In the same section we shall consider the contribution of the spin-factor. As a 
point of note it should be mentioned that, in order to compare our results with standard phe- 
nomenology, we shall adopt some lattice-based parametrizations of the two-point correlator. 
The technicalities of this matter will be discussed in an Appendix. 

2. Semiclassical Calculation 

According to Eq (3), in order to obtain the full amplitude it does not suffice to determine 
the minimal surface bounded by a, given, specific contour -a problem which, in general, is 
very hard to solve 1 . One needs to proceed even further and sum over all possible boundaries 
with a weight of the form 



The particular method we shall follow for getting an estimate of the scattering amplitude 
is the minimization of the correlator contribution to the action, i.e. A[C], related to the 
contour C, while at the same time find the minimal surface corresponding to this specific 
boundary. In this way one obtains, in principle, a result which enables one to determine the 
dominant contribution to the path integral (3). 

1 Indeed, known cases for which solutions have been found are limited and rather simple. 




(6) 



o 
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As was proved in I the variation of A[C] under changes of the boundary reads: 

^ (r) = ^T = ia(T) / dS ^ z ')^W)-z{r\s')} (7) 

^ ' 5(C) 

Accordingly, the correlator contributions become stationary for the "classical" trajectory 

gM = 0. (8) 

Using the expansion for the correlator according to Eq(27) of the previous paper, it is 
easy to see that 

1 

g^T) = 2x a (T)R a ^[x\ - -x a {T)Q aiJi [x\, (9) 

with 

R«M= I dS a ,(z')D[x(r)-z(s',r')) (10) 

S(C) 

and 

QaM = J dr' [x^t') (x a (r') - x a (r)) - (pt <-> a)] D 1 [x(r) - x(r% (11) 

where the quantities D and D x are introduced in the framework of the SVM [1]. Their 
significance is of practical importance, as far as the credibility of the SVM is concerned. Our 
eventual numerical estimates in this work will use them as basic input. It should be further 
noted that the above expressions are reparametrization invariant. Also, in the last relation 
the integration covers the whole range of the r variable. Our next step is to specify the 
minimal surface relevant to the problem under study. 

Following Ref [3] the minimal surface bounded by two infinite rods at a relative angle 9, 
has(in four- dimensional Euclidean space) the shape of a (three-dimensional) helicoid, which 
is the only surface that can be spanned by straight lines. In the considered process the 
eikonal lines 1 — > 3, 2 — > 4, play the role of the 'rods', while the angle 9 is connected, via 
analytic continuation, to the logarithm of the incoming energy. 

Given the above specifications, consider the following, helpful, parametrization of the 
boundary C: For < r < T\ we have a straight line segment, a;' 1 ', going from the point 
Xi to the point X2- Introducing, moreover, for convenience the length 2T =\ x 4 — X2 \ and 
reparametrizing according to r — > — T, we write 

x« = (r, 0,0,0), —T < t <T, (12) 

5 



with xW(-T) = x A) if'(T) = x 2 . 

The second eikonal line x^(r), r 2 < r < r 3 , goes from the point X\ to the point x 3 at 
a relative angle 9 with respect to x^ l \ while a distance 6 (impact parameter) separates the 
two linear contours in a transverse direction. Introducing the distance 2T X =| x 3 — x\ | and 
reparametrizing according to 

r - T X ( -?—t - (13) 



T3 — t~2 r 3 — r 2 



we write 



xj, 3) (r) = (-tcos#, -rsin#, 6, 0), -7\ < r < 71, 



with af)(-Tx) = an, ^(Tx) = x 3 . 

In the following we shall assume, just for convenience, that 

2T =1 X4 — X2 |~ I ^3 — xi 1= 2Ti. 



For T\ < r < r 2 , we have a helical curve xffi(r), which joins the points x 2 = x^(ti) and 
X\ = x^(r 2 ), representing the exchanged light quarks. Performing, now, the change 
s = — - — (t — T\), we write 

4, ( ^(^)cos^,( S )3 i n- s ,0),0< S <, (14) 
The continuity of the boundary requires 

X W(T) = x^(0) = x 2 and = x^(-T) = x u 

or 

0(0) = 0(5) = T. (15) 

The final helical curve is x^(t), which, for r 3 < r < T4, joins the points a; 3 = x^(r 3 ) and 
X4 = x^fa). Making one more, final, parametrization, namely s = (t — r 3 ) we write 



x ji ( s ) = ( -0(s)cosy, -0(s)sin— , s,0 j , < s < b. (16) 



Once again, Eq. (15) takes care of the continuity of the boundary. Now, the minimal surface 
is bounded by the (four) curves specified by Eqs, (12)-(16) and can be spanned by straight 
lines parametrized as follows 

= ^4 4) (*) + ^4 2) (*) = (^0(s)cos^, ^0( S )sin^, s,o) . (17) 



It can be easily proved that the surface defined by the above equation is minimal, irrespec 
tively of the function 0: 

(z ■ z') Zfl - z 2 z'^ 



3 T 



(Z ■ Z')z'^ - Z% 



+ d s 



V9 



0. 



(18) 



One observes that the minimization of the surface is not enough for the complete spec- 
ification of the parametrization of the helicoid. Accordingly, we go back to Eq. (8), which 
determines the boundary that dominates the path integration (3). A first observation is 
that, due to the antisymmetric nature of R afl and Q afM , the function vanishes when x M (r) 
represents a straight line. Thus Eq. (8) is trivially satisfied for the eikonal sector of the 
boundary. Non-trivial contributions are coming only from the helices xff and xff . One can 
simplify Eq. (9) by computing the leading behavior of the functions i? QAt and Q aii using the 
fact that the functions D and Di, as defined in the SVM scheme -and measured in lattice 
calculations [5]- decay exponentially fast for distances which are large in comparison with 
the correlation length T g [1]. In this connection and upon writing 

x(s') = x(s) + (s' - s)x(s) + hs' - s) 2 x(s) H , 



we find, for the second term in Eq. (9), 

1 



[x )x M 



[X ■ X)X, 



J ds'{s' 



X 



7^2 
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+ ■■■ = 



X 



X ■ X 



X — 



ar 



T„ai 



+ 



(19) 



where 2 



oo 

— = Tt\ dww 2 D 1 (w 2 ). 
a± J 



Noting that 



z lt (s,r = T) = x^(s), zJs,r = —T) = xWs) 



(4) | 



fj, v~/j ■ - / ~\i 

1 

2T 1 



d T z^(s,T) = i M (s,r) = — [x {2) (s) -x^\s)], 



(20) 



2 We have omitted terms suppressed by powers of X? 
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the leading behavior of the first term of the rhs of (9) can be easily determined. One finds 



x a R, 



a/i 



1 . 2 / . (x • Z 

2 X { Z> *~^~' J 



T b 

dr' [ ds'D 







x 2 
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+ 



= 2T I x 



{^z)_ u 



x' 



(21) 



where we have introduced the parameter 

oo 

fi 2 = T 2 J dwD{w 2 ). 



Thus, the function g takes, to leading order, the form 



(22) 



1 



57* = 



x T n 



4Tfi 2 x 2 ( z 



x ■ z . 



X ■ X . 



(23) 



X* ' J 2cti \ ' X' 

Now, we recall from its definition, cf Eq. (7) that the (^-function provides a measure of 
the change of A[C] when the Wilson contour is altered as a result of some interaction which 
reshapes its geometrical profile. In this sense, it contains important information concerning 
the dynamics of the problem under study. The structure of the ^-function, as it appears in 
the above equation, is quite general and exhibits its dependence, not only on the boundary 
but on the minimal surface as well. It is worth noting that this fact is strictly associated with 
the non-Abelian nature of the theory since the function D -and consequently fi 2 - disappears 
[1] in an Abelian gauge theory. 

Taking into account that for the helicoids parametrization the velocity x has three non- 
zero components, while 'x and z have only two, we conclude that Eq. (8) can be satisfied 
only if 



2«i 



= 0. 



(24) 



Inserting in Eq. (24) the helical parametrization one easily finds that the function must 
be a constant. Taking, now, into account Eq. (14) we determine this constant to be the 
length T. It is then very easy to see that this result leads to the conclusions 



x ■ z — 0, x ■ x — 



and 



x 



i e 2 

8fi 2 a\ b 2 ' 



(25) 



(26) 



where 

T 2Q2 



x 2 = 1 + 



b 2 

This equation cannot be satisfied in Euclidean space. In Minkowski space the angle 9 
becomes imaginary 9 — > — ix — In (^?) an< ^ Eq. (26) has a positive definite solution: 

TV 1 y 2 

6 2 8/i 2 ai 6 2 ' 1 J 

The above formula relates the impact parameter b, the logarithm of the incoming energy 

In (2^?) anc ^ the distance T. These parameters must not be considered as independent from 

each other in a calculation of the leading behavior of the scattering amplitude. In fact, 

Eq. (27) indicates that the effective impact parameter must grow with the incoming energy: 

b ~ Ins, a conclusion which is in agreement with the landmark result of Cheng and Wu [6]. 

The preceding analysis, obviously repeats itself for the two helical curves x^ and x^ 

and has led us to a specific parametrization for the Wilson loop, which plays the dominant 

role in the path integration in Eq. (3). We are now in position to determine the leading 

contribution to the action (6): 



1 



S cX = - A I drx 2 d {r) + A[C] d . (2£ 



Our first step is to expand the second term of the integrand in powers of T 2 yK. The first 
term of such an expansion is the familiar Nambu-Goto string. The next term, which reveals 
the rich structure of the SVM, is the so-called 'rigidity term', representing the extrinsic 
curvature of a surface embedded in a four-dimensional [7] background: 

A[C] = a f d 2 ^ +— f d 2 z^g ab d a t^d b t^ + ■■■, (29) 

where, in the above expression, 

a^ l -T 2 J d 2 zD(z 2 ) (30) 

enters as the string tension. 

The coefficient of the rigidity term reads 

^ = h T *J d2zz2{2Di{z2) - D{z2)) - (31) 



Terms proportional to T ff 6 entering the expansion in Eq. (28) will be considered negligible 
in our analysis. We have also omitted the term / d 2 £y/gR, since in two dimensions the 
curvature is a total derivative. Using the helicoids parametrization (17), with <fi = T, the 
Nambu-Goto term in Eq. (29) takes the form 

T b I 



J d 2 ^ = JdrJ d S] 



'1 + 



r 2 9 2 



bT 



^l+p 2 + ^ In (^Jl +p 2 + pj 



(32) 



where p — nr- 



To proceed further we analytically continue to Minkowski space where we can use Eq. 
(27) to determine 

i x 2 \ 1/2 1 



bTJl+p 2 -> bT\ 1 - b I 1 - 



b 2 



8a lf i 2 b 2 J v^T? ~ v^T? + ) (33) 



and 

bT_ 

P 

Thus 



ln(vW + p)--^ln 



b 2 1 b 



Tib 2 



2x vSoiJ? 



+ <D(T*). (34) 



° J d 2 ^ 



aub 2 



(35) 



In the same framework, the contribution of the rigidity term takes the form 



/ d 2 CVgg ab d a t^d b t^ 



= 9 



, , , l (e 2 id 4 2 

ar / as , = — H — — r 
1 I U 2 2 6 4 

3, / r \ 1 



-In Ul+p 2 +p) + -p\j 1 + p 



(36) 



It follows that in Minkowski space we have 

l~ J d 2 ^g ah d a t, v d b t 



3vr 



-x- 



(37) 



For the full estimation of the classical action, c/Eq. (28), one should also take into account 
the presence of the classical kinetic term. Non trivial contributions come from the helical 
curves x (2) (l -> 2) and x^(3 -> 4): 



4(r 2 - n) 



I ds(x^f 



+ 



4(r 4 - r 3 ) 



V 

J ds(x^) 2 = 



b 2 x 2 b 2 x 2 

+ 



4(r 2 - ri) 4(r 4 - r 3 ) ' 



(38) 
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Now we have to take into account that both r 2 — T\ and r 4 — r 3 must be integrated with weights 
e -(r 2 -Ti)m 2 anc j e -(r 4 -r 3 )m 2 ^ reS p ec ti V ely. These integrals, as it turns out, are dominated by 
the values r 2 — T\ = r 4 — r 3 = leading to a final kinetic contribution of the form 

2m6 | £ |= 2 v. (39) 

Here, m is the mass of the light quarks, thus the result expressed by (34) can be considered 
negligible. 

From the above analysis we conclude that 

anb 2 37T / s \ 

g ci~ - 7 ; x - — In — - . (40) 



Putting aside, for now, the possible corrections to A[C] which arise from fluctuations of the 
boundary as well as the spin factor contribution, let us consider the result (39) as a whole, 
except for terms ~ m. To obtain the final expression for the scattering amplitude one must 
integrate over the impact parameter: 



/ 



d 2 bexp (iq- b- ^b 2 j oc exp (~^:<? 2 X ) • (41) 



Combining (40) and (41) we find, for the scattering amplitude, a Regge behavior of the 
form s <(o)t +aR (o) with 

1 3n 
°M0) = ^ a*(0) = — . (42) 

In the Appendix we present a certain parametrization [8] for the functions D and Di 
entering the SVM scheme which give for the string tension the value a ps 0.175 GeV 2 and 
for the coefficient of the rigidity term the value ^ 0.276. With these numbers we obtain 
for the Regge slope the value a^(0) 0.91 GeV" 2 and for the Reggeon intercept the value 
aij(0) ~ 0.65, in good agreement with the phenomenological values ot' R (0) = 0.93 GeV and 
a R = 0.55. [3] 

3. Boundary Fluctuations and the Role of the Spin Factor 

As repeatedly mentioned in our narration, corrections to the amplitude (3), beyond 
semiclassical ones, are expected to arise from fluctuations of the boundary of the surface on 
which the two-point correlator 'lives'. Fluctuations of the surface itself can be taken into 
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account by higher order correlators. This, in fact, is the big difference which distinguishes 
the SVM approach, in comparison with Nambu-Goto type approaches. 

We begin our related considerations by expanding the action (6) around the helicoid 
classical solution: 

S = S d - - J dTy(T)x cl (r) + - j dr J dfy a r) x 



x 





d 2 







1 . d 2 .. 5 2 A[C] 
—2 5a ^ 5{T - T) + 5xi{r)5x {f) 



(43) 



where y = x — x d . 

Using the results of I one can easily determine that 

5 2 A[C] 



5x a (T)5xp(f) 



x,{r)x v {f)A^[x{r)-x{f)] - 



J dS Xp (z')AW Xp [z(e-x(T)} + 



S(C) 



+x a {r) J dsa{f,s)z x {f,s)z' p {f,s)e KuXp A Kafl [z{f,s) -x{t)], (44) 



where we have written 



5x v 



= ShJ(t - f)a(f,s). 



The second term on the rhs of Eq.(44) is simply the area derivative which, as we have seen 
in I, has the general form ~ g a £p — 9px a - Thus, for the classical solution g[x cl ] it 

gives zero contribution. It is, furthermore, easy to verify that the third term in (44) also 
disappears for x = x cl . We, therefore, conclude that 

5 2 A[C] 



x 



■^T)x*{f)A% p [x c \r)-x{f) 



(45) 



5xi{r)5x${f) 

Inserting Eq. (45) into Eq. (44) and taking into account that the dominant contribution 
to the two-point correlator comes from the region r pa f we find 



S pa S cl + J dsy a (s) 



1 m d 2 X 2 



2 x 



ds 2 T n 



(46) 



Let it be remarked that to arrive at the above relation we have adopted the expansion of 
the two-point correlator indicated in Eq. (21) of paper I. We have also used the helicoid 
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parametrization observing, at the same time, that the eikonal lines give null contribution. 
One further realizes that the contributions of the two helical curves to the linear term in 
(42) cancel each other, since x^ 2 \s) = —x^(s) and r-i — T\ ~ T4 — T3 ~ 

The non-trivial contribution of the helical curves is incorporated in the term 

^ = S aP -^(x?xf+x^xf), (47) 

the origin of which is the second functional derivative, c.f. (48). The mass parameter A 2 in 
(46) has the same source and is defined as 



X 2 =\x\ T 2 J dw (d(w 2 ) + D 1 (w 2 ) + A.D 1 (w 2 )) . (48) 

The differential operator entering Eq. (46) has no zero eigenvalues since the "classical" 
solution is, in fact, the one that annihilates the (^-function. Accordingly, the calculation 
of the path integral over y = x — x cl does not require any particular regularization. A 
straightforward calculation shows that 

If 1 \ T 2 6 2 /b 2 . . 

det ^ = ^( 1 -^) = iT^W- (49) 

Thus the matrix oj a p can be diagonalized and the y-integral can be easily performed. How- 
ever, in the limit m — > it can be immediately seen that the integration over the boundary 
fluctuations gives prefactors which are powers of the logarithm of the incoming energy and 
as far as Regge behavior is concerned, they cannot change the behavior that was determined 
in the previous section. 

The next task is to take up the issue of the spin-field dynamics contribution to the 
scattering amplitude. As seen in I, a spin factor is associated with each segment of the 
worldline path. This factor receives contributions from two sources. The first one is 

Jdr J dS- A( 2 \z -x).J = J dr ^-^ i [ 7/l , 7 „] (50) 

5(C) 

and is obviously zero for the classical trajectory (8) 
The other term has the form 

1 /■,/■.,,. o\ , 3 



S =8 



Jdr J dr'J^A {2 l Xp (x - x')J Xp = - J ' dr J dr'(D + D ± ) + - J dr J dr'(x - x') 2 D[. 

(51) 
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In the stochastic limit, within which we are working, the integrals in the above equation 



give appreciable contribution to (51) only for | x(r) 



X(T 



X 



t-t' |< T g . More con- 



cretely, consider the contribution to (51) from the helical curve (1 — > 2). A straightforward 
calculation shows that the analytically continued result is 



S=-(t 2 - tl f— , 

X 

where we have written t — it for the time variable and denoted 

\ 1/2 

on 



(52) 



M 



\ o 



/ dwD(w) 
_o 

/ dw w 2 D 1 (w) 



1 d 



dw\D{w 2 )+D l {w 2 ) + -—D l {w 



(53) 



As has been mentioned in I and discussed in [9], contribution (52) has an interesting role 
as far as the form of the fermionic propagator is concerned, but it is obvious that it does not 
alter the basic Regge structure of the amplitude was calculated in the previous section. 
The remaining spin structure is summarized in the chain 

1 



rrii 



i=4 



2m- 



(54) 



which must be sandwiched between the external spinor wavefunctions representing the in- 
coming and outgoing quarks (in the simple picture wherein the meson wavefunction is just 
the product of free spinors). The non-trivial dynamics of the process are now incorporated 
into the fact that the vectors xff, % — 1, 2, 3, 4 forming the boundary of the helicoids, are 3- 
dimensional vectors with | | 2 = const. For i — 1,3 turns the factor in (54) to the operator 



X2 



(i) 



\ p (i)\ ■ 

For i = 2,4 the matrices 
/ 2 = 1 



and 



L = 1 



2m r 2 — T\ 
1 b 



1 -x i2 \b) -> 1 



7-^ 4 )(6) 



1 - 



1 -x {2 \b) 

I i( 2 ) I 



7 -±( 4 )(6) 
I W) I : 



(55) 



(56) 



2m T4 — r 3 

are also representations of projection operators. As shown in [3] the matrices (55) and (56) 
are the direct product of two 2x2 matrices each of which are by themselves projection 
operators. Given these observations it becomes a matter of simple algebra to find that the 
standard kinematics are reproduced. 
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Appendix 



Appendix 

In this Appendix we present a parametrization of the functions D and Di, already referred 
to in I and used extensively in the present paper. This parametrization is supported by 
lattice data and is extensively discussed in Ref. [8]. 
The exact relations defining the functions are 

2N C 24 ' 2N C 24 v ; ' ' v ; 

where D N and D 1N are functions which determine the structure of the two-point correlators, 
as defined in [8]. The factor G 2 is defined as follows 

„2 r) AT 

G 2 ee (0 | |-i£(0)F« (0) | 0) = ^A^(0). (A.2) 

For the above correlator we shall adopt the value given in Ref [8], namely G 2 = (0A96) A GeV 4 . 
The value of the numerical quantity k in (A.l) is estimated in the same reference to be 0.74 

The ansatz for the function D N is [8] 

D ^> = I i / d,ke ""'- r ~- — JT4 > (A.3) 

where 

oo 

a = J dzD N (z). (A.4) 
o 

A simple calculation shows that 

1 3n 
D N (z) = wKi(w) - -w 2 K (w), w = — \ z \, (A. 5) 

4 oa 
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with K v denoting a Bessel function. The correlation length T g , introduced in I -and fre- 
quently used in the text- can be deduced from Eq. (A. 5): 

T, = -. (A.6) 

The estimated value of a is 

a « 0.35 fmor T g « 0.297fm. (A.7) 
With the help of ansatz (A. 3) and using (A.7) one can determine the string tension: 



V f d 2 wD(w) = -T f {N j- l) % I 



wK\{w) - -w 2 K (w) 



or 



(A.8) 



2 9 J v ' 2 9 2N C 24 

or 

a = \ T9 ^n~ 1) Ta k27{ = a2G ^fi * °- 175Gey2 (A ' 9) 

The anzatz for the function -D^jvis deduced from the equation [8] 



U + ZlM -^j D hN (z) = AD N (z) (A.10) 



D liN (z) = \\ Z dw\Aw A K x {w) - w 5 K (w)] (A.ll) 
z Jo 

The coefficient of the rigidity term entering Eq. (31) can now be calculated: 
^ = ^Jd 2 ww 2 [2D 1 (w)-D(w)] = 

= h ^wc 1 ^ I rfW[2(1 - K)DiAw) - kDn{w)] 

= ^ T f {I li~ l) ^{l - n)32n » 0.276. (A.12) 
32 9 2N C 24 v 1 v 1 
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